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Abstract. We describe the implementation of a biochemi- 
cal model of isoprene emission that depends on the elec- 
tron requirement for isoprene synthesis into the Farquhar- 
Ball-Berry leaf model of photosynthesis and stomatal con- 
ductance that is embedded within a global chemistry-climate 
simulation framework. The isoprene production is calculated 
as a function of electron transport-limited photosynthesis, 
intercellular and atmospheric carbon dioxide concentration, 
and canopy temperature. The vegetation biophysics mod- 
ule computes the photosynthetic uptake of carbon dioxide 
coupled with the transpiration of water vapor and the iso- 
prene emission rate at the 30 min physical integration time 


step of the global chemistry-climate model. In the model, 
the rate of carbon assimilation provides the dominant con- 
trol on isoprene emission variability over canopy tempera- 
ture. A control simulation representative of the present-day 
climatic state that uses 8 plant functional types (PFTs), pre- 
scribed phenology and generic PFT-specific isoprene emis- 
sion potentials (fraction of electrons available for isoprene 
synthesis) reproduces 50% of the variability across differ- 
ent ecosystems and seasons in a global database of 28 mea- 
sured campaign-average fluxes. Compared to time-varying 
isoprene flux measurements at 9 select sites, the model au- 
thentically captures the observed variability in the 30 min 
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average diurnal cycle ( R 2 = 64-96 %) and simulates the flux 
magnitude to within a factor of 2. The control run yields 
a global isoprene source strength of 451TgCyr _1 that in- 
creases by 30 % in the artificial absence of plant water stress 
and by 55 % for potential natural vegetation. 


1 Introduction 

Global terrestrial gross primary productivity (GPP) is the to- 
tal amount of carbon dioxide (CO 2 ) removed from the atmo- 
sphere by plant photosynthesis. GPP is the largest flux in the 
global carbon cycle and is of fundamental importance to life 
on Earth as the basis for food and fiber. By absorbing an in- 
creasing amount of fossil fuel CO 2 , GPP provides a critical 
ecosystem service of climate protection in the Anthropocene 
(Ballantyne et al., 2012). Estimates of GPP center around 
120PgCyr _1 (Beer et al., 2010) although recent isotopic 
analysis supports a higher value of 150-175 PgCyr -1 (Welp 
et al., 2011). Land ecosystems return to the atmosphere an 
estimated 1 % of GPP in the form of biogenic volatile or- 
ganic compounds (BVOCs) (Guenther et al., 2012). This 
> 1 PgC yr ~ 1 chemically reactive carbon flux is of compa- 
rable magnitude to the annual global net ecosystem produc- 
tion and an order of magnitude larger than the annual anthro- 
pogenic VOC source. The ecological and physiological roles 
of BVOCs are broad and range from abiotic and biotic stress 
functions to integrated components of carbon metabolism 
(Loreto and Schnitzler, 2010; Kesselmeier and Staudt, 1999). 
The dominant BVOC emitted is isoprene, amounting to half 
of the total annual flux of reactive carbon (Guenther et al., 
2012). The biogenic isoprene flux is comparable to the total 
flux of methane including both biogenic and anthropogenic 
sources. The specific roles of isoprene emission are not fully 
understood but may be related to protection against heat 
and antioxidants. For instance, isoprene stabilizes chloro- 
plastic membranes during high temperature events, allowing 
the plant’s photosynthetic capacity to be maintained during 
rapid leaf temperature fluctuations caused by sun flecks in the 
canopy (Sharkey and Singsaas, 1995; Behnke et al., 2010). 
Isoprene reduces plant damage caused by ozone and reactive 
oxygen species (Vickers et al., 2009). 

The chemical transformation of isoprene in the atmo- 
sphere has a profound effect on the distribution and vari- 
ability of the key short-lived climate forcers: tropospheric 
ozone (O 3 ), secondary organic aerosol, and methane (CH 4 ) 
(Fiore et al., 2012). Thus, the isoprene flux is of central im- 
portance to understanding interactions between atmospheric 
chemistry and climate. Indeed, isoprene emission may have 
played an important role in the Earth system climate sen- 
sitivity in past greenhouse worlds (Beerling et al., 2011). 
While it is qualitatively perceived that human and/or natu- 
ral perturbations to isoprene emission may provide a power- 
ful lever on regional climate and even trigger feedbacks to 


global climate (Pitman et al., 2012), there is a lack of quan- 
titative understanding of this critical ecosystem-chemistry- 
climate linkage. The main challenges have been the accu- 
rate model representation of the isoprene emission response 
to the complex, sometimes, opposing, influences of global 
change (Monson et al., 2007) and the development of ap- 
propriate coupled global carbon-chemistry-climate modeling 
frameworks. 

Isoprene emission rate depends strongly upon ecosystem 
type with broadleaf trees and shrubs exhibiting the strongest 
emission potentials. Isoprene is produced in the chloroplast 
from precursors formed during photosynthesis. Isotopic la- 
beling shows that about 70-90% of isoprene production is 
directly linked to photosynthesis, which provides the sup- 
ply of energy, reducing power and carbon skeletons for in- 
leaf biosynthesis (Delwiche and Sharkey, 1993; Karl et al., 
2002; Affek and Yakir, 2003). The remaining isoprene pro- 
duction is associated with an older carbon source. Any sin- 
gle molecule of isoprene can include precursors from both 
newly assimilated and the older carbon source. In the labora- 
tory and the field, the incomplete coupling to the photosyn- 
thetic flow is revealed in 4 ways: (1) the time lag between 
the onset of photosynthesis and isoprene emission, which 
is due to the effects of the growth environment on the ex- 
pression of isoprene synthase e.g. (Kuzma and Fall, 1993; 
Sharkey and Loreto, 1993; Fuentes et al., 1999; Monson et 
al., 1994; Goldstein et al., 1998; Pressley et al., 2005); (2) the 
higher temperature optimum of isoprene production versus 
photosynthesis such that at high leaf temperatures the rates 
of isoprene production and photosynthesis are inversely cor- 
related; (3) the initial short-term increase in isoprene emis- 
sion under drought stress e.g. (Pegoraro et al., 2004, 2005, 
2006). Under sustained drought stress, the isoprene emission 
rate does decrease according to the reduction in photosyn- 
thetic rate. Active upregulation of isoprene production dur- 
ing drought stress has been posited (Monson et al., 2007; 
Niinemets, 2010) but there remain unresolved issues regard- 
ing the isoprene emission response to drought (4) the iso- 
prene emission rate has an inverse relationship in response 
to increasing atmospheric CO 2 concentration (opposite to 
that observed for the photosynthetic CO 2 assimilation rate), 
behavior known as the “CO 2 -inhibition effect” (Rosenstiel 
et al., 2003; Possell et al., 2004, 2005; Possell and Hewitt, 
2011; Monson et al., 2007). Competition for carbon sub- 
strate has been offered as a mechanistic explanation for the 
CO 2 -inhibition effect (Rosenstiel et al., 2003; Wilkinson et 
al., 2009). Enhanced isoprene-related tolerance of heat and 
light stressed photosynthesis has been found at low but not 
high CO 2 concentrations (Way et al., 2011). It has been hy- 
pothesized that isoprene biosynthesis evolved in a low CO 2 
climate state, which prompts the important question: what is 
the functional benefit of leaf isoprene production in a future 
high CO 2 world? 

Global isoprene emission models have been developed 
for application in chemistry-climate studies (Guenther et al., 
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1995, 2006; Lathiere et al., 2006; Arneth et al., 2007; Paci- 
fico et al., 201 1). The most widely used approach to simulate 
interactive isoprene emission is to modify ecosystem-specific 
basal emission rates under standard conditions, either at the 
leaf or canopy level, with empirical functions in the form 
of serial multipliers that describe the observed emission re- 
sponse to specific environmental controls (light, temperature, 
soil moisture) (Guenther et al., 1991, 1995, 2006). This ap- 
proach has been applied within prescribed static vegetation 
frameworks, for instance in the complex global canopy en- 
vironment Model of Emissions of Gases and Aerosols from 
Nature (MEGAN) (Guenther et al., 2006); and within dy- 
namic global vegetation model (DGVM) frameworks, for in- 
stance in the Organizing Carbon and Hydrology in Dynamic 
EcosystEms model (ORCHIDEE) (Lathiere et al., 2006) and 
in the Community Land Model (Heald et al., 2008). Sim- 
ilarly, the empirical functions that describe isoprene emis- 
sion response to light and temperature at the leaf or canopy 
level have been directly embedded within global chemistry- 
climate and global chemistry-transport models (CCMs and 
CTMs) using prescribed static vegetation e.g. (Wu et al., 
2008; Shindell et al., 2006; Horowitz et al., 2007). Current 
generation global CCMs and CTMs usually neglect the re- 
sponse of isoprene emission to soil moisture. Precipitation 
controls photosynthesis in more than 40 % of vegetated land 
(Beer et al., 2010). It follows that water availability must play 
an important role in isoprene emission that will in turn affect 
the atmospheric composition and climate-air pollution inter- 
actions. 

Medium to long-term changes of vegetation physiology 
and composition in response to global change drivers (cli- 
mate change, CO 2 , land use), and increasing awareness of 
the complex mutual feedbacks between isoprene emission 
and regional climate sensitivity, make it necessary to link 
isoprene emission directly to the biological processes that 
affect emissions. To address this urgent need, a leaf-level 
isoprene emission model that depends on the electron re- 
quirement for isoprene synthesis (Niinemets et al., 1999) 
has been modified for implementation within DGVM frame- 
works including the Lund Potsdam Jena General Ecosys- 
tem Simulator (LPJ-GUESS) (Arneth et al., 2007), and the 
Joint UK Land Environmental Simulator (JULES) (Pacifico 
et al., 2011). Isoprene emissions generated in this way have 
been used in chemistry-climate modeling studies. For in- 
stance LPJ-GUESS isoprene emissions were applied off-line 
in the UM CAM model to quantify the effects of the CO 2 - 
inhibition on future ozone predictions (Young et al., 2009), 
and JULES isoprene emissions have been applied on-line in 
the HadGEM2 Earth-system model to examine the sensitiv- 
ity of isoprene emission to past and future global change and 
the implications for atmospheric chemistry (Pacifico et al., 
2012 ). 

Global modeling of carbon cycle-climate and chemistry- 
climate and interactions have evolved as entirely sepa- 
rate communities because CO 2 is chemically unreactive 


in the atmosphere and because of the vast differences in 
the system integration time scales. There are many ex- 
isting DGVMs whose contemporary global carbon cycle 
simulations have been extensively evaluated (Sitch et al., 
2008; Schwalm et al., 2010). These models do not typi- 
cally include on-line atmospheric chemistry. On the other 
hand, global chemistry-climate models normally rely on 
off-line prescribed vegetation input data sets. Our objec- 
tive is to implement the biochemical leaf isoprene produc- 
tion scheme of (Niinemets et al., 1999) into on-line vegeta- 
tion biophysics integrated within a global chemistry-climate 
model framework. To achieve this objective, we introduce 
the Yale-E2 global carbon-chemistry-climate model that is 
built around the NASA Goddard Institute for Space Stud- 
ies (GISS) Model-E2 global climate model and features ad- 
vanced biogeochemistry-climate interactions under separate 
development at Yale University. 

The major goals of this study are (i) to describe the model 
and (ii) to evaluate the global scale model performance in the 
present climate state. The modeling methodology and iso- 
prene emission algorithm are detailed in Sect. 2. In Sect. 3, 
we specify the control and sensitivity global simulations that 
are carried out for this work. Section 4 summarizes the global 
model results and presents the model evaluation including: 
summary of simulated global GPP and isoprene emission 
(Sect. 4.1), isoprene emission sensitivity to GPP and canopy 
temperature (Sect. 4.1.1), evaluation of global GPP and iso- 
prene emission (Sect. 4.2), FLUXNET-derived global GPP 
(Sect. 4.2.1), GPP and latent heat seasonal cycle at 6 bench- 
mark sites (Sect. 4.2.2), global database of campaign-average 
isoprene flux measurements (Sect. 4.2.3), time-varying iso- 
prene emission through campaign periods (Sect. 4.2.4), iso- 
prene and GPP diurnal cycle (Sect. 4.2.5). Discussion and 
conclusions are presented in Sect. 5. 


2 Methodology 

2.1 Yale-E2 global carbon-chemistry-climate model 

Yale-E2 is built around the new generation IPCC AR5 ver- 
sion NASA Model-E2 global climate model (Schmidt et al., 
2006) and incorporates interactive terrestrial ecosystems, a 
dynamic carbon cycle module, and 2-way coupling between 
the on-line vegetation and atmospheric chemistry. The model 
has flexible horizontal and vertical resolution. In this study, 
we apply 2° x 2.5° latitude by longitude horizontal resolu- 
tion with 40 vertical layers extending to 0.1 hPa. The vege- 
tation submodel is embedded within the general circulation 
model that provides the key meteorological drivers for the 
vegetation physiology (Friend and Kiang, 2005). The land- 
surface hydrology submodel provides the grid cell level soil 
characteristics to the vegetation physiology. The well estab- 
lished gas-phase chemistry and aerosol modules are fully in- 
tegrated, so that these components interact with each other 
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Table 1. PFT-specific fractional coverage of global vegetated land area in the standard and SiB2 data sets (%) and isoprene emission 
parameters used in Yale-E2 global carbon-chemistry-climate model. PFT-specific photosynthesis parameters are in Table Al. 


PFT 

SimCONT fraction (%) 

SimSIB2 fraction (%) 

Mw gCg *h ') 

£ 

Tundra 

5.5 

8.1 

0 

0.000 

Grass 

6.4 

5.4 

16 

0.016 

Shrub 

3.2 

8.2 

16 

0.055 

Savanna 

4.6 

16.7 

16 

0.036 

Deciduous 

10.5 

7.4 

45 

0.061 

Evergreen 

7.8 

7.8 

8 

0.015 

Rainforest 

9.8 

9.7 

24 

0.027 

Crop 

36.9 

19.1 

0 

0.000 


and with the physics of the climate model (Bell et al., 2005; 
Shindell et al., 2006, 2013; Unger, 201 1). 

2.1.1 Vegetation structure 

The vegetation is described using 8 plant functional types 
(PFTs): tundra, grass, shrub, savanna, deciduous, tropical 
rainforest, evergreen, and crop (Table 1). In this work, we 
apply two different vegetation cover data sets that have been 
converted to the 8 PFTs: (i) the standard atlas-based distribu- 
tion in NASA Model-E2 (Matthews, 1983) and (ii) the Sim- 
ple Biosphere Model II (SiB2) distribution based on the In- 
ternational Satellite Land Surface Climatology Project data 
initiative II (Loveland, 2009). The SiB2 data set provides its 
own crop cover. The standard data set consists of a map of the 
world’s vegetation cover that would most likely exist in equi- 
librium with present-day climate and natural disturbance, in 
the absence of human activities (potential natural vegetation) 
onto which the crop fraction for each grid cell is overlaid. 
The crop fraction in each model grid cell is constructed from 
a harmonized gridded data set for the year 2000 (Hurtt et 
al., 2011). The crop cover is imposed by proportional de- 
crease/increase of all the potential natural vegetation types in 
the grid cell fraction that is not occupied by crops and/or pas- 
ture. There is no right choice in how to implement the crop 
cover, but this approach is the most common treatment of 
crop cover in global climate modeling (de Noblet-Ducoudre 
et al., 2012). The fractional coverage of global vegetated land 
area by each PFT for the standard and SiB2 data sets is shown 
in Table 1 (bright and dark bare soil fractions not shown). 
In the standard data set, the fractional cover of crop PFT is 
about double that in SiB2 (36.9% versus 19.1 %) because 
tropical and subtropical land that is classified as crop PFT in 
the standard data set is classified as savanna and shrub PFTs 
in SiB2. 

Leaf area index (LAI) for each PFT is prescribed accord- 
ing to regular seasonal sinusoidal variation between PFT- 
specific minimum and maximum seasonal LAI values that 
is insensitive to climate drivers or carbon balances (Rosen- 
zweig and Abramopoulos, 1997; Friend and Kiang, 2005). 
A complete mechanistic understanding of the processes that 


control the development and senescence of foliage is not yet 
available such that the current state of phenology modeling 
may even be considered qualitative (Migliavacca et al., 2012; 
Richardson et al., 2013). We have made some improvements 
to the phenology for this work. Firstly, we have implemented 
a simplified crop phenology using a global data set of crop 
planting and harvesting dates (Sacks et al., 2010). The dom- 
inant crop type in each model grid cell was identified using 
a published data set of global crop maps and areal cover- 
age (Monfreda et al., 2008). Then, a global model input file 
of mean plant date and harvest date was constructed for the 
grid cell’s dominant crop type. The plant and harvest dates 
are recycled every simulation year. Secondly, a parameteri- 
zation for phenological control (frost hardening) on photo- 
synthetic capacity (Vcmax, Sect. 2.1.1) has been added for the 
evergreen PFT. Vcmax is reduced in winter to protect against 
cold injury. Sensitivity to temperature determines the revival 
of photosynthetic capacity, hence the length of the growing 
season (Toivonen et al., 1991; Makela et al., 2004; Hanninen 
and Kramer, 2007). The onset of the growing season for the 
evergreen PFT can be delayed if temperature remains cold 
despite light being available. 

The use of fixed canopy structures and phenology means 
that leaf mass is not driven by photosynthetic uptake of CO 2 
and a closed carbon cycle is not simulated. However, this ver- 
sion of the isoprene emission model can respond to elevated 
CO 2 with regard to: CO 2 fertilization, reduced stomatal con- 
ductance/increased water-use efficiency, CO 2 -inhibition, the 
temperature and precipitation responses of photosynthesis, 
but not phenological timing because the simulation does not 
account for the coupling of photosynthetic uptake to vari- 
ability in growth. Application of LAI that is insensitive to 
climate may dampen the simulated interannual variability of 
isoprene emission in this model version. 

2.1.2 Canopy biophysics 

Each model PFT fraction in the vegetated part of each grid 
cell represents a single canopy. The model vertically strati- 
fies each canopy into diffuse and direct light levels, and LAI 
profiles using an adaptive number of layers (typically 2-16) 
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(Friend and Kiang, 2005). The well established Michealis- 
Menten leaf model of photosynthesis (Farquhar et al., 1980; 
von Caemmerer and Farquhar, 1981) and the stomatal con- 
ductance model of Ball and Berry (Collatz et al., 1991) 
is used to compute the biophysical fluxes at the leaf level 
in each canopy layer based on appropriate parameters for 
each of the 8 PFTs from (Friend and Kiang, 2005) and the 
Community Land Model (Oleson et al., 2010) with updates 
from (Bonan et al., 2011) (Table Al). This coupled pho- 
tosynthesis/stomatal conductance leaf model has previously 
been widely used to project terrestrial biosphere responses 
to global change. We summarize the model briefly here for 
transparency and completeness. The photosynthesis model 
assumes that the rate of net CO 2 assimilation in the leaves 
of C3 plants is limited by one of three processes: (i) the ca- 
pacity of the ribulose 1,5-bisphosphate (RuBP) carboxylase- 
oxygenase enzyme (Rubisco) to consume RuBP (Rubisco- 
limited photosynthesis); (ii) the capacity of the Calvin cycle 
and the thylakoid reactions to regenerate RuBP supported 
by electron transport (electron transport-limited photosyn- 
thesis); and (iii) the capacity of starch and sucrose synthe- 
sis to consume triose phosphates and regenerate inorganic 
phosphate for photo-phosphorylation (triose phosphate use- 
limited photosynthesis). Photosynthesis is electron transport- 
limited under low light conditions including overcast/cloudy 
conditions, at the start and end of the day, for shaded leaves 
and understory vegetation. The three processes are described 
as functions of the internal leaf CO 2 concentration (C/) 
and/or the maximum carboxylation capacity at the optimal 
temperature, 25 °C, (Vcmax) (Table Al). Leaf stomata control 
for the uptake of CO 2 versus the loss of water vapor (H 2 O). 
In the model, the stomatal conductance of H 2 O through the 
leaf cuticle is linearly related to the net rate of carbon assim- 
ilation and the relative humidity, and inversely related to the 
CO 2 concentration at the leaf surface. The coupled system of 
photosynthesis, stomatal conductance and CO 2 diffusive flux 
transport equations form a cubic in C; that is solved analyti- 
cally (Baldocchi, 1994). A simple but realistic representation 
of soil water stress is included in the vegetation biophysics 
following the approach of Porporato et al. (2001). The algo- 
rithm reflects the relationship between soil water amount and 
the extent of stomatal closure ranging from no water stress to 
the soil moisture stress onset point (s*) through to the wilt- 
ing point (s w iit). Stomatal conductance is reduced linearly be- 
tween the PFT-specific values of s* and s w ii t (Table 1) based 
on the climate model’s soil water volumetric saturation in 6 
soil layers. 

The leaf-level carbon and water fluxes are scaled up to 
the canopy level by integrating over each canopy layer. The 
land-surface model uses its own internal adaptive time step 
from 5 s to 15 min depending on the conditions. The carbon 
and water fluxes computed in the land-surface scheme are 
integrated over the climate model time step (30 min) for ex- 
change with the model’s atmosphere. 


2.1.3 Leaf isoprene production 


A leaf-level isoprene emission model that describes the con- 
stitutive production as a function of the electron transport- 
limited photosynthesis rate, J Q , (Niinemets et al., 1999) has 
been integrated into the canopy biophysics scheme follow- 
ing modifications for global-scale modeling (Arneth et al., 
2007). The leaf-level isoprene emission rate (I) in units of 
pmol m -2 [leaf] s -1 is calculated as follows: 

/ = J e • ft • K • T • £, (1) 


where J Q is the electron transport limited photosynthesis rate 
in units of jumol m“ 2 [leaf] s“ 1 . J Q is a linear function of the 
incident photosynthetically active radiation (PAR) and the in- 
ternal leaf CO 2 concentration (Ci): 


J Q — d\ Qa f • PAR • Qfqe • 


et - r* 
ci - 2r* ’ 


( 2 ) 


where a\ ea f is the leaf-specific light absorbance and a qQ is the 
intrinsic quantum efficiency for photosynthetic CO 2 uptake 
in the chlorophyll reaction system that absorbs PAR to drive 
the oxidation of water and the reduction of enzymes (photo- 
system II). a qQ is a product of the fraction of absorbed light 
that reaches photosystem II and the CO 2 per absorbed pho- 
ton. T* is the CO 2 concentration compensation point in the 
absence of non-photorespiratory respiration (Collatz et al., 
1991). 

The p term in Eq. (1) translates the electron flux into iso- 
prene equivalents given by Eq. (3): 


P = 


Ci - r* 

6(4. 67C; + 9.33T*) 


(3) 


A detailed description of the mechanistic origin of the co- 
efficient values is given elsewhere (Niinemets et al., 1999; 
Pacifico et al., 2011). 

The atmospheric CO 2 -inhibition is included via a simple 
parameterization (k): 


Ci standard 



(4) 


where Ci standard is the leaf internal CO2 concentration at 
standard atmospheric CO 2 , which is chosen to be the year 
2000 global average value (370ppmv). Equation (4) mim- 
ics the observed response to both short-term and long-term 
changes in C/ (Wilkinson et al., 2009; Heald et al., 2009). 
For example, short-term reductions in Q due to stomatal clo- 
sure under drought conditions imply increases in k, while 
the long-term effects of increasing (decreasing) atmospheric 
CO 2 imply decreases (increases) in k. 

The temperature relationship (r) in the algorithm accounts 
for the difference in temperature optimum between photo- 
synthesis and isoprene synthase: 


r = exp[o.i(r-r ref )], 


(5) 
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where T = leaf temperature and r re f= standard temperature 
condition (30 °C). The temperature optimum for isoprene 
synthase is about 40 °C (Guenther et al., 1991). Equation (5) 
does not simulate a temperature optimum after which iso- 
prene emission rate decreases with further increases in tem- 
perature. Such high temperature conditions in isoprene emit- 
ting biomes rarely occur in nature at large ecosystem scales. 
Canopy-scale temperatures of this magnitude may occur un- 
der severe drought stress conditions when transpiration is 
significantly reduced. Yale-E2 uses the canopy temperature 
(not air temperature) in Eq. (5), which represents a significant 
departure from other CCMs and CTMs that drive interactive 
isoprene emissions with surface air temperature. Yale-E2 in- 
trinsically captures the effects of changing stomatal conduc- 
tance on canopy energy balance, which affects the canopy 
temperature, and thus the isoprene emission rate. In future 
and past hot greenhouse worlds, plant photosynthesis may 
acclimate to the higher temperatures (Ameth et al., 2012). 
This plastic adjustment of photosynthesis will indirectly im- 
pact isoprene emission. Whether the temperature optimum 
for isoprene synthase will similarly shift in warmer climates 
is not known. 

s is the PFT-specific fraction of electrons available for iso- 
prene synthesis or isoprene emission potential and parallels 
the use of a PFT-specific leaf-level isoprene emission capac- 
ity in Guenther et al. (1991, 1995). To calculate e, we make 
use of available generic PFT-specific standard (or basal) leaf 
isoprene emission rates (7 S ) in jugC g — 1 [leaf] h -1 based on 
recommendations from a wide range of observations (Ta- 
ble 1). Standard conditions are defined as: surface air tem- 
perature = 30°C and PAR = 1000 pmolm -2 s -1 . Again, we 
assume that atmospheric CO 2 concentration = 370 ppm un- 
der standard conditions. 7 S must be converted to appropri- 
ate units of leaf area using the model values of specific leaf 
area (SLA), the amount of light-capturing surface area that 
is deployed with a given investment of dry mass in units of 
m 2 [leaf] g _1 [leaf]. SLA is intimately connected to the re- 
source use economy of the plant (Milla and Reich, 2007). 
Then, e is computed for each PFT by substituting 7 = 7 S into 
Eq. 1) under standard conditions. The crop and tundra PFTs 
are non-emitting for isoprene in this model although recently, 
isoprene emission has been observed from a tundra ecosys- 
tem (Potosnak et al., 2013). Deciduous, shrub and rainforest 
PFTs maintain the highest fraction of electrons available for 
isoprene synthesis (Table 1). In the current model, s does 
not vary with time of day or season (Niinemets et al., 2010a, 
c) or through the canopy (Niinemets et al., 2010b). Many 
plant species in temperate and boreal ecosystems exhibit de- 
layed onset of isoprene emission after the leaf development 
(Grinspoon et al., 1991; Hakola et al., 1998; Olofsson et al., 
2005). Leaf age effects have been included in other global 
isoprene emission models using simplistic functions, for ex- 
ample involving growing degree-day temperature (Ameth et 
al., 2007; Guenther et al., 2006; Lathiere et al., 2006). 


Scaling of isoprene emission from the leaf to the canopy 
uses the canopy vertical stratification and integration scheme 
as described in Sect. 2.1.2. The canopy-level isoprene fluxes 
are passed to the model’s atmosphere through the land- 
surface scheme on the 30 min climate model time step. 

3 Simulations 

We apply Yale-E2 in an atmosphere-only configuration 
driven by sea surface temperatures and sea ice cover pre- 
scribed according to decadal average monthly varying obser- 
vations for 1996-2005 from the HadSST2 data set (Rayner et 
al., 2006). Atmospheric CO 2 is prescribed to a uniform con- 
centration of 370 ppmv representative of approximately year 
2000 levels. We perform a control simulation (“SimCONT”) 
that uses the standard atlas-based vegetation cover distribu- 
tion in NASA Model-E2 (Matthews, 1983). In addition, we 
perform three sensitivity simulations described in Table 2. 
Sensitivity simulation “SimSiB2” uses an alternate vegeta- 
tion cover data set based on the SiB2 distribution. We per- 
form a simulation “SimH20” in which the plant water stress 
function in the biophysics module is artificially switched off 
to quantify the effects of water limitation on isoprene emis- 
sion rate. In order to explore the extent of human land cover 
change impacts on isoprene, a simulation “SimPNV” is run 
that uses the potential natural vegetation cover data set (i.e. 
without the crop PFT superimposed). The control and sen- 
sitivity simulations are each run for 10 model years that are 
averaged for the evaluation analyses. The model version ap- 
plied in this study outputs the total isoprene (and GPP) in 
each grid cell and does not characterize PFT-specific attri- 
bution. The model’s carbon fluxes (GPP and isoprene) are 
present-day climatologies and do not pertain to any specific 
meteorological year. 

4 Results 

4.1 Simulated global GPP and isoprene emission 

The global annual average isoprene emission and GPP flux 
for each simulation are reported in Table 2. Global GPP in 
SimCONT and SimSiB2 (124-126 PgCyr -1 ) is in reason- 
able agreement with current understanding of the contempo- 
rary carbon cycle budget. The global isoprene source in Sim- 
CONT and SimSiB2 (451-498 TgCyr -1 ) is consistent with 
previous global estimates of 400-700 TgCyr -1 (Guenther et 
al., 2006). This good agreement is because the model ap- 
plies widely used generic PFT-specific 7 S values to derive the 
PFT-specific fraction of electrons for isoprene synthesis ( s ) 
as previously alluded (Ameth et al., 2008). The spatial dis- 
tribution of annual and seasonal average GPP and isoprene 
emission are shown in Fig. la, b, respectively. The major 
isoprene emitting biome is tropical rainforest that provides 
a perennial reactive carbon flux to the atmosphere, but the 
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Table 2. Summary of global annual average carbon fluxes in the model runs. In total 10 model run years are included in the annual average 
for each simulation. All runs use 1996-2005 average monthly varying SSTs and sea ice. Uncertainty values are based on the standard error 
relative to internal climate model variability ( n = lOyr). 



Simulation Land cover 


Description and purpose 


GPP 


Isoprene 


(PgCyr l ) (TgCyr l ) 


SimCONT 

Standard, Matthews (1983) 

Present Climate Control 

124 ±2 

45 1 =b 5 

SimSiB2 

SiB2, Loveland (2009) 

Test impact of alternate vegetation cover data set 

126 ±2 

498 ±5 

SimH20 

Standard, Matthews (1983) 

Test impact of water stress function 
in vegetation biophysics 

155 ±1 

592 ±3 

SimPNV 

Standard + no human alteration, Matthews (1983) 

Test impact of potential natural vegetation 

137 ± 2 

701 ±10 


ANNUAL 


ANNUAL 



Fig. la. Spatial distribution of annual and seasonal average modeled Fig. lb. Spatial distribution of annual and seasonal average modeled 

GPP in gC m“ 2 day -1 in the control run SimCONT. isoprene emission in mgC m“ 2 h" 1 in the control run SimCONT. 
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growing season isoprene emission in the eastern US rivals 
the tropical emission in magnitude. 

Globally, the isoprene source strength is 10% higher in 
SimCONT than SimSiB2 whereas GPP is only ~1 % higher. 
The different response sensitivity of isoprene and GPP to the 
vegetation cover data set arises because of the differences in 
the classification of crop PFT between data sets, and all PFTs 
assimilate carbon but the crop PFT is non-emitting for iso- 
prene. Previous studies have suggested an even higher sensi- 
tivity of the global isoprene source (20-30 %) to the model’s 
vegetation cover data sets for a given algorithm (Guenther et 
al., 2006; Pfister et al., 2008; Ameth et al., 2011). In the ar- 
tificial absence of plant water stress, GPP increases by 25 % 
with a similar increase in global isoprene emission (^30 %). 
Similarly, the MEGAN model suggested that accounting for 
soil moisture stress in the present climate state leads to a 
20-30% reduction in the global isoprene source strength 
(Muller et al., 2008). The response to potential natural veg- 
etation cover is more dramatic: the global isoprene source 
increases by 55% whereas GPP increases by only 10%. 
The strong sensitivity of isoprene to potential natural veg- 
etation cover suggests significantly higher isoprene emission 
in pre-human disturbance conditions and large variability in 
isoprene emission due to human land cover change. In the 
most recent community assessment of anthropogenic ozone 
radiative forcing led by the Atmospheric Chemistry and Cli- 
mate Model Intercomparison Project, only 4 of 15 partici- 
pating state-of-the-science global chemistry-climate models 
included climate-sensitive isoprene emission (Young et al., 
2013). Those 4 models all projected a small increase in iso- 
prene emission from preindustrial to present-day in response 
to temperature change. 9 out of 15 models prescribed iso- 
prene emission and used the same off-line input data for 
preindustrial and present-day. 

There are substantial regional differences in isoprene 
emission between the sensitivity and control simulations 
(Fig. 2). Annual average isoprene emission in Europe and 
the southeastern US is smaller in SimCONT than SimSiB2 
by ~0.5— 1 .5 mgC m -2 h _1 . In the tropical biome, SimCONT 
simulates higher annual average isoprene emission than Sim- 
SiB2 by ~0.3-l mgCm -2 h -1 . The presence of plant water 
stress affects isoprene in all emitting biomes. In the artifi- 
cial absence of plant water stress, large increases in the an- 
nual average isoprene emission (~1 mgCm -2 h -1 ) occur in 
the tropics, southeastern US, and the boreal region. Potential 
natural vegetation cover implies additional isoprene emission 
(~2 mgCm -2 h -1 on the annual average) in the eastern US, 
Europe and the tropics on land that has zero emission capac- 
ity in the standard cover data set. 

4.1.1 Isoprene emission sensitivity to GPP and canopy 
temperature 

We have applied multiple linear regression to compute the 
standardized regression coefficients between isoprene emis- 



- 2.1 -1 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1 1.7 

SimH 20 - SimCONT 




- 1.2 -1 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1 2 


Fig. 2. Difference in simulated annual average isoprene emis- 
sion relative to the control simulation for an alternative vegeta- 
tion cover map (top), in the artificial absence of water stress (mid- 
dle), and for potential natural vegetation cover (bottom). Units 
are mgCm -2 h -1 . 


sion and the key drivers (GPP and canopy temperature) 
in boreal summer for SimCONT and SimE120 (Fig. 3). 
The coefficients represent the partial regression coefficients 
in units of standard deviation, such that they can be di- 
rectly compared with each other to determine the rela- 
tive importance of the different driving variables (Fig. 3). 
The leaf-level isoprene emission rate is linearly related 
to the electron transport-limited photosynthesis rate, J Q 
(Eq. 1). Above 30 °C, the isoprene emission follows an 
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GPP 


Canopy T 








- 0.9 - 0.7 - 0.5 - 0.3 - 0.1 0.1 0.3 0.5 0.7 0.9 1.1 1.3 1.5 


Fig. 3. Standardized regression coefficients between isoprene emission and the key drivers GPP (left column) and canopy temperature (right 
column). Results shown are averages for June- July- August in SimCONT (top row), SimH20 (middle row), and the difference [SimH20 - 
SimCONT] (bottom row). 



Latitude (°) 


Fig. 4. Seasonal zonal average GPP simulated in SimCONT 
and SimSiB2 compared to the FLUXNET-derived data set 
(gCm -2 day -1 ). 


exponential relationship with temperature (Eq. 5). In the 
control run (SimCONT), the standardized regression coef- 
ficient is about a factor of 2 higher for GPP than canopy tem- 


perature everywhere except the most severely water-limited 
ecosystems (e.g. Indian subcontinent). Thus, GPP dominates 
over canopy temperature in controlling isoprene emission 
variability in most regions of the world in the model. In the 
absence of water stress (SimH20), the standardized regres- 
sion coefficients for canopy temperature increase by around 
30 % while the standardized regression coefficients for GPP 
decrease by around 30 %, relative to the control run (Sim- 
CONT). However, GPP still exerts the dominant control over 
canopy temperature in controlling isoprene emission vari- 
ability in this model. The presence of plant water stress re- 
duces the sensitivity of isoprene emission to canopy temper- 
ature because drought conditions are frequently associated 
with high temperatures that result in a decrease in GPP and 
therefore isoprene emission. Two important corollaries are 
that (1) variability in isoprene emission is largely determined 
by the same factors that determine GPP variability (tempera- 
ture, precipitation, the availability of water and nutrients, dis- 
turbance) and (2) isoprene emission algorithms that do not 
account for the impacts of water availability may overesti- 
mate the sensitivity of isoprene emission to canopy (and/or 
surface air) temperature. Since a global observational data 
set of isoprene emission is not available, we are unable to 
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explore these relationships at large regional scales in the real 
world. 

4.2 Evaluation of global GPP and isoprene emission 

Direct large-scale measurements of photosynthesis and iso- 
prene emission do not exist. The FLUXNET network of 
tower sites continuously measures the CO 2 net ecosystem 
exchange (NEE) flux between the biosphere and atmosphere 
using the eddy covariance technique (Baldocchi et al., 2001). 
NEE is the imbalance between photosynthesis and ecosystem 
respiration. Photosynthesis is zero at night, which allows for 
a temperature-based predictor of the respiration that can be 
used to extrapolate GPP from the day time NEE. Isoprene 
fluxes have been measured during short campaigns at several 
FLUXNET sites. Our approach to evaluate the global model 
performance is to compare the simulated GPP and isoprene 
emission rate to available above-canopy flux tower measure- 
ments of the carbon fluxes. Recent research has shown that 
space-based observations of formaldehyde (HCHO) columns 
may be a useful proxy of surface isoprene emissions across 
broader ecosystem and regional scales (Palmer et al., 2003; 
Barkley et al., 2009). We elect not to use this approach here 
because of the inherent limitation that an a priori isoprene 
emission model is necessary to interpret the HCHO satellite 
columns (Barkley et al., 2012). 

Because of the different spatial and temporal scales of the 
model output versus the observations, comparing to point 
measurements from flux towers represents an extremely 
stringent performance test and a number of caveats must 
be emphasized. Firstly, the model GPP and isoprene emis- 
sion are extracted for the single grid cell (2° x 2.5° lati- 
tude by longitude) in which the flux tower site is located. 
Our rationale is that, while the typical eddy covariance foot- 
print is small (~1 km 2 ), the trace gas flux variability at the 
flux towers is often representative over much larger spatial 
scales because of the spatial coherence of climate anoma- 
lies (Ciais et al., 2005). Secondly, the model parameters are 
not tuned in any way to the local site vegetation properties. 
The model grid cell level output is based entirely on the 
vegetation structure and soil properties in Yale-E2 and not 
the flux tower site data. Finally, the point measurements are 
obtained at flux towers during specific meteorological peri- 
ods. The model output represents lOyr average climatolog- 
ical GPP and isoprene emission rates (and meteorological 
variables) at the site locations, which do not necessarily re- 
flect the exact local weather conditions that occurred during 
the observation period. Since the climate model has not been 
forced to the observed meteorology during the measurement 
periods, the model cannot be expected to reproduce day-to- 
day variability in the isoprene emission and GPP. Specifi- 
cally for isoprene, the measurements are canopy exchange 
fluxes whereas the current model configuration does not in- 
clude isoprene loss through the canopy. However, canopy 
loss is likely less than 10% of the total emission (Karl et 


al., 2004; Rinne et al., 2012). Based on the above caveats, 
we posit that the most important diagnostic quantities in this 
model/measurement comparison are: zonal average and sea- 
sonal cycle for GPP flux, and campaign-average variability 
across different ecosystems and diurnal cycle for isoprene 
emission. Focusing on these diagnostic quantities does pro- 
vide for valuable insights into the global climate model’s 
strengths and weaknesses. Typical measurement uncertain- 
ties associated with the eddy flux technique are about ±30 %. 

4.2.1 FLUXNET-derived global GPP 

We compare the model simulated GPP (SimCONT and Sim- 
SiB2) to a global GPP data set that has been generated 
by data orientated diagnostic upscaling of site-derived GPP 
from FLUXNET (Beer et al., 2010; Jung et al., 2011; Bonan 
et al., 2011). The model runs reproduce the seasonal zonal 
average variability in the FLUXNET-derived data set with 
remarkable fidelity (Fig. 4). Based on linear regression, the 
model runs capture > 89 % of the variability in the zonal sea- 
sonal average FLUXNET-derived data set. The model runs 
perform best in the summer and winter seasons (> 93 % of 
the variability). SimSiB2 demonstrates slightly superior per- 
formance over SimCONT. The Northern Hemisphere (NH) 
June- July- August maximum at 50°-60°N in the decidu- 
ous and crop biomes is well reproduced by the model (to 
within 10%). In March- April-May and September-October- 
November, the model underestimates tropical GPP and over- 
estimates GPP in the NH mid-latitudes between 35°-60°N 
by around 10 %. Overall, the comparison results provide con- 
fidence in the model’s ability to simulate the magnitude and 
zonal average variability in GPP. 

4.2.2 GPP and latent heat seasonal cycle at 
6 benchmark sites 

A benchmarking system of seasonal FLUXNET data at se- 
lect sites has been constructed for evaluation of the seasonal 
cycle of carbon and water biophysical fluxes in global mod- 
eling frameworks (Blyth et al., 2011). The characteristics of 
the 6 benchmark sites that are relevant for isoprene emis- 
sion are described in Table 3. At the Hyytiala, Harvard and 
Santarem sites, the dominant vegetation PFT in the model 
grid cell is consistent with the FLUXNET tower site. At Mor- 
gan Monroe and Tharandt, the dominant model PFT in the 
grid cell is crop versus deciduous and evergreen, respectively 
at the FLUXNET tower sites. The Kaamanen FLUXNET 
tower site is wetland/woody savanna in reality whereas the 
standard cover is predominantly grass PFT in this model grid 
cell, and the SiB2 cover is a mixture of evergreen, deciduous 
and shrub PFTs. 

The climatological models capture the GPP seasonality at 
the 6 sites reasonably well (Fig. 5a). We provide root mean 
square bias (RMSE) as the diagnostic metric (Table 3). All 
model climatologies overpredict at both NH high latitude 
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Table 3. Summary of FLUXNET benchmark site characteristics, model vegetation cover fractions for the grid cell in which the site is located 
and model performance in terms of RMSE range across the model runs. Td = tundra; G = grass; Sh = shrub; S = savanna; D = deciduous; 
T = tropical rainforest; E = evergreen; and C = crop. 


Name 

Data year 

Site characteristics 
Ecosystem Type 

Location 

Vegetation cover fraction (%) 
SimCONT SimSiB2 

GPP 

(jimol m 2 s" 

RMSE 

Latent Heat 
_1 ) (Wm“ 2 ) 

Kaamanen 

2002 

Wetland/woody savanna 

69.1° N, 27.3° E 

G:D 80/10 

E:Sh:D:Td:S:C 40/31/17/9/3 

1.3-1. 5 

7.4-9. 0 

Hyytiala 

2000 

Evergreen needleleaf 

61.8° N, 24.3° E 

E:D:C 35/34/31 

E:C:D 75/19/5 

1.4— 1.7 

12.7-14.4 

Morgan Monroe 

2002 

Deciduous broadleaf 

39.3° N, 86.4° W 

C:D 51/49 

C:D 67/31 

1.7-1. 8 

20.5-21.6 

Harvard 

1999 

Deciduous broadleaf 

42.5° N, 72.2° W 

D:C 97/3 

D:E:C 79/16/5 

1. 9-2.6 

24.0-25.6 

Tharandt 

1999 

Evergreen needleleaf 

51.0° N, 13.6° E 

C:E:D 51/25/24 

C:D:E 82/15/2 

2.2-3A 

14.4-16.4 

Santarem Km67 

2003 

Tropical evergreen broadleaf 

2.9° S, 55.0° W 

T:C 94/6 

T:S:E 87/7/3 

0.4— 1.0 

6.1-10.5 



Month Month 


Fig. 5a. Comparison of eddy-derived and modeled monthly average 
GPP fluxes at 6 FLUXNET sites (pmolC m -2 s _1 ). 


sites (Kaamanen and Hyytiala) and underpredict at NH mid- 
latitude sites (Morgan Monroe, Harvard and Tharandt). The 
model demonstrates skillful ability to reproduce the GPP flux 
at the tropical site (Santarem). SimCONT and SimSiB2 cap- 
ture the observed GPP reduction in the dry season that is over 
predicted in the artificial absence of water stress (SimH20). 
The model’s ability to reproduce the magnitude of summer 
GPP at Harvard is improved in the absence of water stress 
(SimH20) but that configuration gives too high GPP in the 
spring and fall. 

Previous analyses of FLUXNET measurements at a wide 
range of sites found a strong linkage between carbon gain and 
water loss for all biomes integrated over the year (Law et al., 
2002). Figure 5b compares the modeled and measured latent 
heat (evaporative) fluxes at each of the 6 FLUXNET sites. 
RMSE values are shown in Table 4. Yale-E2 satisfactorily 



Month Month 


Fig. 5b. Comparison of measured and modeled latent heat flux at 6 
FLUXNET sites (Wm“ 2 ). 

captures the seasonal cycle in latent heat flux at all sites, ex- 
cept at Harvard. Again, artificially turning off the water stress 
function does appear to improve the model’s latent heat sim- 
ulation at the Harvard site. The model’s high values of latent 
heat outside of the growing season are due to the grid cell 
having 25 % water coverage. 

The model/measurement comparison of the GPP seasonal 
cycle at the 6 benchmark sites permits reasonable confi- 
dence in the model’s biophysics for the simulation of iso- 
prene emission. Based on the model GPP biases alone, the 
model may be expected to overestimate isoprene emission at 
NH high latitudes, underestimate at NH mid-latitudes, and 
perform well in the tropics. 


www.atmos-chem-phys.net/13/10243/2013/ 


Atmos. Chem. Phys., 13, 10243-10269, 2013 


10254 


N. Unger et al.: Photosynthesis-dependent isoprene emission from leaf to planet 


Table 4. Description of the 28 campaign-average above-canopy isoprene flux measurements collected between 1995-2010 across a wide 
range of ecosystem types, regions and seasons that comprise the benchmark global database. Codes for model vegetation fractions are: 
Td = tundra; G = grass; Sh = shrub; S = savanna; D = deciduous; T = tropical rainforest; E = evergreen; and C = crop. 


Identifier 

Ecosystem 

Measurement sites 

Date 

Location 

Reference 

Model grid cell vegetation fractions 
Standard SiB2 

Santarem 

Rainforest 

07/2000 

2.9° S, 55° W 

Rinne et al. (2002) 

T: 94 

T/Sv/E: 87/8/4 

Manaus 

Rainforest 

07/2001 

2.35° S, 60° W 

Kuhn et al. (2007) 

T: 99 

T/Sv: 89/5 

Manaus 

Rainforest 

09/2004 

2.6° S, 60° W 

Karl et al. (2007) 

T: 99 

T/Sv: 89/5 

Costa Rica 

Rainforest 

04/2003 

10° N, 84° W 

Karl et al. (2004) 

T/C: 82/18 

T/C: 74/14 

Costa Rica 

Tropical lowland 

10/1999 

10° N, 84° W 

Geron et al. (2002) 

T/C: 82/18 

T/C: 74/14 

Malaysia 

Rainforest 

11/2003 

3° N, 102° E 

Saito et al. (2008) 

T/C: 79/21 

C/T/D: 45/27/16 

Borneo 

Rainforest 

04/2008 

5° N, 118° E 

Langford et al. (2010a) 

T/C: 87/13 

T/E: 78/15 

Borneo 

Rainforest 

06-07/2008 

5° N, 118° E 

Langford et al. (2010a) 

T/C: 87/13 

T/E: 78/15 

Congo 

Mixed tropical 

03/1996 

2° N, 16° E 

Serca et al. (2001) 

T/C: 93/7 

T/95 

Congo 

Mixed tropical 

11/1996 

2° N, 16° E 

Serca et al. (2001) 

T/C: 93/7 

T/95 

Congo 

Mixed tropical 

11/1996 

4° N, 18° E 

Greenberg et al. (1999) 

T/Sv/C: 72/18/10 

T/Sv: 76/19 

South Africa 

Savanna 

02/2001 

25° S, 31° E 

Harley et al. (2003) 

C/Sh/Sv: 64/22/14 

C/D/T: 66/14/18 

Mongolia 

Grassland 

06-07-08/2002 

44° N, 116° E 

Bai et al. (2006) 

C/Sh: 97/2 

G/C: 84/10 

Mongolia 

Grassland 

09/2002 

44° N, 116° E 

Bai et al. (2006) 

C/Sh: 97/2 

G/C: 84/10 

SW China 

Plantation 

02-03/2002 

21° N, 101° E 

Baker et al. (2005) 

T: 100 

T/D/E: 43/27/27 

SW China 

Plantation 

07/2002 

21° N, 101° E 

Baker et al. (2005) 

T: 100 

T/D/E: 43/27/27 

Germany 

Deciduous 

07/2001-2002 

50° N, 11° E 

Graus et al. (2006) 

C/D: 48/47 

C/D: 84/15 

NW Germany 

Deciduous 

07/2003 

50° N, 6° E 

Spirig et al. (2005) 

C/D: 52/48 

C/D: 85/14 

Germany 

Deciduous 

08/1995 

50° N, 8° E 

Steinbrecher et al. (2000) 

D/C: 59/41 

C/D: 74/24 

Sweden 

Plantation 

07/2001 

58° N, 12° E 

Olofsson et al. (2005) 

E/C/D: 49/35/16 

C: 97 

Finland 

Boreal 

08/2001 

61° N, 24° E 

Spirig et al. (2004) 

E/D/C: 35/33/32 

E/C: 75/19 

N Wisconsin 

Deciduous 

07/1993 

46° N, 90° W 

Isebrands et al. (1999) 

D/C/Sv: 43/37/19 

C: 86 

UMBS 

Deciduous 

06-07-08/2000-2002 

46° N, 85° W 

Pressley et al. (2005) 

D/E/C: 57/29/13 

E/D/C: 36/33/31 

UMBS 

Deciduous 

08/1998 

46N, 85W 

Westberg et al. (2001) 

D/E/C: 57/29/13 

E/D/C: 36/33/31 

Harvard Forest 

Deciduous 

06-07-08/1995 

45° N, 72° W 

Goldstein et al. (1998) 

D/C: 92/8 

D/C: 64/26 

Canada 

Boreal 

07/1994 

42° N, 72° W 

Pattey et al. (1999) 

C/E/Sv: 67/16/10 

C/E: 70/13 

Canada 

Boreal 

09/1994 

54° N, 105° W 

Pattey et al. (1999) 

C/E/Sv: 67/16/10 

C/E: 70/13 

London 

Urban 

10/2006 

54° N, 105° W 

Langford et al. (2010b) 

D/C: 57/43 

C: 98 


4.2.3 Global database of campaign-average isoprene 
flux measurements 

Our goal is to assess the model’s ability to capture the iso- 
prene emission variability across a wide range of ecosys- 
tem types and seasons. We have assembled a global database 
of campaign-average above-canopy isoprene flux measure- 
ments obtained between 1995-2010 (Table 4). The measure- 
ments are sorted by ecosystem type. The dominant vegeta- 
tion fraction in the model grid cells for both the standard 
and SiB2 vegetation cover data sets corresponds to the lo- 
cal ecosystem at the measurement towers for all sites in the 
database except for South Africa and SW China. Table 5 
compares the model simulated isoprene emission in Sim- 
CONT, SimH20 and SimSiB2 to the observations in the 
benchmark database (shown graphically in Fig. 6 for Sim- 
CONT and SimH20 results). The model results are monthly 
or seasonal averages that correspond to the time period dur- 
ing which the observations were made. The campaign peri- 
ods are typically shorter than one month and as such the ob- 
servations tend to refer to shorter averaging periods than in 
the model output. In some cases the published observational 
campaign-average is for mid-day average or daytime aver- 


age as indicated in Table 5. Evaluation against this global 
database suggests possible systematic biases in the model: 
under prediction in the summer time Canadian boreal but 
over prediction in the Scandinavian boreal. The model per- 
forms well in simulating the magnitude of isoprene mission 
in the tropical biome but appears to simulate either a lack of 
seasonality or reverse seasonality compared to that observed 
at the non- Amazon tropical rainforest sites. The ORCHIDEE 
model isoprene emission has been compared to a limited sub- 
set of this global database (Canada, UMBS and Santarem) 
(Lathiere et al., 2006). Similar to Yale-E2, ORCHIDEE re- 
produced the tropical flux magnitude well but underpredicted 
at the Canada and UMBS sites. The presence of water stress 
has the largest impact on the magnitude of isoprene emission 
rate in the North American deciduous sites. Based on linear 
regression, the control run SimCONT is able to reproduce 
50 % of the variability across different ecosystems and sea- 
sons in the global database. SimH20 and SimSiB2 capture 
64 % and 40 % of the variability in the global database, re- 
spectively. The model/measurement discrepancies in the iso- 
prene emission magnitude are most likely due to the model’s 
assignment of PFT-specific emission potentials (fraction of 
electrons available for isoprene synthesis, s) that in reality 
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Table 5. Comparison of the global database and simulated isoprene emission. Units are mgC m~ 2 h _ 1 . The model values are either monthly 
or seasonal means to be consistent with the observation time period. The averaging period for the observations is indicated where the 
information is available. [MDA = mid-day average; DTA = daytime average; CA = campaign average; DTMX = daytime maximum]. 


Location 

Measurement 


Model 


Santarem 

1-2 MDA 

SimCONT 
2.02 ±0.09 

SimH20 
1.85 ±0.07 

SimSiB2 
1.96 ±0.09 

Manaus 

2.1 DTA 

1.00 ±0.42 

1.98 ±0.09 

0.53 ±0.48 

Manaus 

7.8 MDA 

2.58 ± 0.10 

2.51 ±0.09 

2.62 ±0.12 

Costa Rica 

1.35 MDA 

2.09 ±0.22 

2.22 ±0.07 

1.30 ±0.32 

Costa Rica 

2.5 

1.48 ±0.11 

1.67 ±0.12 

1.63 ±0.12 

Malaysia 

1.20 CA 

1.62 ±0.08 

1.59 ± 0.12 

1.22 ±0.05 

Borneo 

0.54 CA 

2.23 ±0.08 

2.06 ±0.09 

1.99 ± 0.10 

Borneo 

1.2 CA 

1.42 ±0.26 

2.00 ±0.10 

1.05 ±0.31 

Congo 

1.4 MDA 

2.02 ±0.10 

1.89 ±0.09 

2.08 ±0.08 

Congo 

0.46 MDA 

1.80 ±0.06 

1.93 ± 0.13 

1.72 ±0.09 

Congo 

0.6-1. 3 

1.33 ± 0.14 

1.78 ±0.09 

1.70 ± 0.12 

South Africa 

0.34 MDA 

0.72 ±0.12 

0.92 ±0.08 

0.49 ±0.21 

Mongolia 

0.5 

0.03 ±0.00 

0.04 ±0.00 

0.48 ±0.05 

Mongolia 

0.38 

0.02 ±0.00 

0.02 ±0.00 

0.30 ±0.02 

SW China 

0.15 DTA 

0.12 ± 0.08 

0.66 ±0.11 

0.06 ±0.15 

SW China 

1 DTA 

0.24 ±0.54 

1.64 ± 0.16 

0.34 ±0.39 

Germany 

0.21-0.33 DTMX 

1.37 ± 0.33 

1.52 ± 0.15 

0.49 ±0.12 

NW Germany 

3.4 CA 

1.36 ± 0.12 

1.68 ± 0.15 

0.48 ±0.05 

Germany 

0.8 CA 

0.95 ±0.26 

1.34 ± 0.17 

0.12 ±0.08 

Sweden 

0.44 CA 

1 .79 ±0.11 

2.26 ±0.19 

0.43 ± 0.04 

Finland 

0.035 CA 

0.78 ± 0.10 

0.97 ±0.08 

0.40 ±0.08 

N. Wisconsin 

1-1.89 DTA 

1.75 ±0.50 

2.20 ±0.15 

0.35 ±0.08 

UMBS, US 

3 MDA 

0.35 ±0.11 

2.56 ±0.14 

0.21 ±0.05 

UMBS, US 

9-10 (warm period) 

0.40 ±0.13 

2.35 ± 0.18 

0.24 ±0.06 

Harvard 

0.9-1. 8 (cool period) 
1.16 DTA 

0.79 ±0.14 

3.40 ±0.14 

0.66 ±0.24 

Canada 

2.29 

0.27 ±0.09 

0.31 ±0.01 

0.40 ±0.10 

Canada 

0.54 

0.13 ±0.03 

0.12 ±0.02 

0.16 ±0.02 

London 

0.13 CA 

0.35 ±0.05 

0.42 ±0.05 

0.01 ±0.00 


vary greatly for a single plant species (Goldstein et al., 1998) 
and across plant species lumped into each PFT class e.g. 
(Rinne et al., 2009). For example, the deciduous PFT in- 
cludes species that are strong isoprene emitters (willows, as- 
pen) and others that do not emit isoprene (birches). 

4.2.4 Time-varying isoprene emission through 
campaign periods 

We examine the time-varying performance at 9 select mea- 
surement sites described in Table 6. Some of the sites are part 
of the global database of campaign-average fluxes (Table 3). 
There are 4 tropical sites on 2 different continents, 4 tem- 
perate broadleaf sites on 2 different continents, and 1 tem- 
perate mixed site in Europe. At 3 of the temperate sites, mea- 
surements are available across multiple years: Harvard Forest 
(1995 and 2007); UMBS (2000, 2001, 2002, 2003, 2005); 
and Belgium (2009-2011). These three sites provide mea- 
surements that span an entire growing season. The remain- 
ing 6 sites offer measurements over periods of a few days 


to several weeks. The sites in southern France (La Verdiere 
and Montmeyan) occur in the same model grid cell. We focus 
on analysis of the control run SimCONT and the sensitivity 
simulation SimH20. Sub-daily model output is generated at 
half-hourly resolution (the physical time step of the climate 
model itself). The model isoprene emissions are lOyr clima- 
tological averages for present-day, thus cannot be expected 
to reproduce the day-to-day weather-related variability in the 
observations. 

In this section, we assess qualitatively the 
model/measurement comparison of isoprene emission, 
surface air temperature (SAT) and downward shortwave 
radiation (SW) across the entire campaign period for the 
6 sites with measurement periods shorter than a growing 
season. The goal is to identify possible model biases in 
isoprene emission magnitude before a more quantitative 
assessment of the average diurnal cycle is executed in 
the next section. Results for tropical sites are shown in 
Fig. 7a-f and results for the temperate sites are shown in 
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Fig. 6. Scatter plot of the simulated isoprene emissions against mea- 
surements from the global above-canopy flux database sorted by 
ecosystem type (Table 5). “Plantation” measurements are lumped 
into “crop” here. Units are mgC m -2 h -1 . Values and temporal av- 
eraging periods are detailed in Table 5. 



Fig. 7a. (a-f) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the tropical measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SimCONT (red) and SimH20 (blue), (a) 
Costa Rica during 17 April-4 May, 2003. 


Fig. 8a-e. Results for Harvard Forest (2007) have been 
split into three separate time periods. Meteorological data 
are not available for all sites. Measured isoprene fluxes are 
highly variable during the day at both tropical and temperate 
sites reflecting the local weather conditions. The model 



Fig. 7b. (a-f) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the tropical measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SimCONT (red) and SimH20 (blue), (b) 
Manaus during 17-27 September 2004. 



Fig. 7c. (a-f) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the tropical measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SimCONT (red) and SimH20 (blue), (c) 
wet season Santarem during 17-21 April 2001. 

reproduces the SW with striking fidelity at both the tropical 
sites and temperate sites (where data is available). The SAT 
simulation is generally weaker and the discrepancies are 
not related to a particular biome. In Costa Rica, the model 
underestimates the diurnal range in temperature. In Manaus, 
La Verdiere and Montmeyan, the model captures accurately 
the daytime maximum SAT but is not able to simulate the 
nighttime minimum (the nighttime SAT is too warm in the 
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Julian Day 


Fig. 7d. (a-f) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the tropical measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SimCONT (red) and SimH20 (blue), (d) 
dry season Santarem during 22 October-6 November 2003. 



Fig. 7e. (a-f) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the tropical measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SimCONT (red) and SimH20 (blue), (e) 
wet season Borneo during 21 April-9 May 2008. 


model versus the observations). In Borneo the model’s SAT 
diurnal range appears reasonable but there is a consistent 
^5 °C positive bias. Only at the Harvard Forest site does the 
model accurately reproduce the observed SAT throughout 
the measurement campaign. It must be emphasized that the 



Fig. 7f. (a-f) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the tropical measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SimCONT (red) and SimH20 (blue).(f) 
early dry season Borneo during 22 June-22 July 2008. 


model’s climatological temperature and radiation variables 
refer to average values over the 2° x 2.5° grid cell. SAT is 
expected to exhibit much higher sub-grid scale variability 
than downward SW. 

The model demonstrates significant skill in closely repro- 
ducing the time-varying isoprene emission in the dry season 
tropics at the Manaus, Santarem and Borneo (early-dry) sites. 
However, there is an obvious model overestimate at the Costa 
Rica, Santarem and Borneo sites in the wet season. A large 
part of the discrepancy is likely related to variation in the 
emission potential (s) across different tropical forests. Some 
of the discrepancy at the Santarem site in the wet season may 
be related to the lack of leaf age effects in the model (i.e. no 
seasonal variation in s). For instance, low HCHO columns 
over Amazonia observed from space during the wet season 
have been ascribed to leaf flushing prior to the dry season 
(Barkley et al., 2009). 

As expected, at the wet season tropical sites, SimCONT 
and SimH20 yield similar climatological fluxes. At the dry 
season tropical sites, especially Santarem and Borneo, Sim- 
CONT reduces to about half the value in SimH20 and more 
closely matches the observations lending some confidence 
in the model’s simulation of water-stressed photosynthesis 
and isoprene emission in this region, previously indicated 
in Fig. 5a. Some of the model overestimate in Borneo may 
be related to the climate model’s high bias in SAT. At the 
temperate sites in the growing season, the isoprene emis- 
sion magnitude is well captured in the model (within a factor 
of 2 of the observations). At Montmeyan (June) both model 
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Fig. 8a. (a-e) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the temperate measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SIMCONT (red) and SimH20 (blue), (a) 
La Verdiere during 21 June-6 July 2000. 


Fig. 8c. (a-e) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the temperate measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SIMCONT (red) and SimH20 (blue), 
(c) Harvard Forest during 13-22 June 2007. 




Fig. 8b. (a-e) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the temperate measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SIMCONT (red) and SimH20 (blue), (b) 
Montmeyan during 1 1-27 June 2001. 


Fig. 8d. (a-e) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the temperate measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SIMCONT (red) and SimH20 (blue), (d) 
Harvard Forest during 20-3 1 July 2007. 


runs perform well, whereas at La Verdiere (June-July), Sim- 
CONT reproduces the measured isoprene flux more closely 
than SimH20, indicative of increases in water stress later 
into the growing season. At the Harvard Forest site, the mea- 
sured fluxes tend to lie between SimCONT and SimH20. 


The model reproduces the decrease in emission magnitude 
in early September at this site. 
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Fig. 8e. (a-e) Time evolution of isoprene emission flux, surface 
air temperature (SAT) and downward shortwave (SW) radiation at 
the temperate measurement sites. Measurements are shown in black 
(where data is available). Model climatological results are averages 
of 10 simulation years in SIMCONT (red) and SimH20 (blue), (e) 
Harvard Forest during 31 August-6 September 2007. 

4.2.5 Isoprene and GPP diurnal cycle 

We constructed campaign average diurnal cycles at 30 min 
resolution for isoprene emission at the 9 measurement sites 
(1 h resolution for Harvard Forest in 1995) and GPP at 2 tem- 
perate sites where suitable data is available (UMBS and Har- 
vard Forest). The measurement days included in the averag- 
ing are indicated in Table 6. At UMBS and Belgium, the iso- 
prene (and GPP at UMBS) average diurnal cycle includes the 
data for all available years. The interannual variability in the 
measured isoprene diurnal cycle (30-100%) is larger than 
for the GPP diurnal cycle (< 5 %) (not shown). 

Comparison between measured and modeled isoprene av- 
erage diurnal cycles at the tropical and temperate sites are 
shown in Fig. 9a and b, and GPP average diurnal cycle at 
2 temperate sites in Fig. 9c. Model/measurement evaluation 
statistics including the correlation coefficient for linear re- 
gression and root mean square error (RMSE) are indicated 
in Tables 7 and 8 along with a comparison of the measured 
and modeled maximum average diurnal emission. The mea- 
sured diurnal average profiles show large fluctuations during 
the middle of the day that are not seen in the smoother 10- 
year average model climatologies. The 5 yr average diurnal 
average at UMBS does not exhibit fluctuations to the same 
extent as the other sites. 

The model simulates the magnitude of the average diurnal 
cycle in the dry season tropics and the growing season tem- 
perate zone (except at Belgium) to within a factor of 2, which 
is comparable to the performance of higher spatial resolution 
models driven with observed meteorology that are used in air 



Fig. 9a. (a-c) Comparison between measured and modeled isoprene 
emission and GPP average diurnal cycles for the time periods in- 
dicated in Table 6. Measurements are black lines. Error bars on 
the measurements represent =b 1 standard deviation. Model results 
are lOyr climatological averages for SimCONT (red) and SimH20 
(blue). The dotted lines represent model uncertainty due to internal 
variability in the climate model computed as =b 1 standard deviation 
for n = 10 model run years, (a) isoprene emission tropical sites (wet 
season results are shown in the left column and dry season results 
are shown in the right column) 


quality modeling for the development of environmental pol- 
icy e.g. (Warneke et al., 2010; Ameth et al., 2011). Based 
on the RMSE (Table 7), the control run SimCONT performs 
better than SimH20 in simulating the emission magnitude, 
especially at UMBS, Harvard Forest and the seasonally dry 
tropics. The large model/measurement discrepancy (^factor 
of 5) at Belgium is because the model classifies the grid cell 
as deciduous PFT whereas in reality the measurement site is 
mixed temperate composed of coniferous (Douglas fir) and 
deciduous (beech). The fraction of electrons available for iso- 
prene synthesis (e) is a factor of 4 higher for the deciduous 
PFT than the evergreen PFT (Table 1). In reality, the contri- 
bution to the isoprene flux at the Belgium site from conifer- 
ous versus deciduous is small but beech is non-emitting for 
isoprene. We retain the Belgium site in this study because it 
illustrates the dependency of the model’s simulation on the 
vegetation data set and structure (Schurgers et al., 2011; Ar- 
neth et al., 2011). 

The simulated isoprene average diurnal cycle is lower at all 
sites (except Manaus) in SimCONT than SimH20 following 
the expected GPP response to the presence of water stress. 
In dry season Manaus, isoprene emission is higher in the 
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Harvard Forest 1995 (June-July-August) 


Harvard Forest 2007 (June-July) 




Montmeyan (June) 


La Verdiere (June-July) 





Fig. 9b. (a-c) Comparison between measured and modeled isoprene 
emission and GPP average diurnal cycles for the time periods in- 
dicated in Table 6. Measurements are black lines. Error bars on 
the measurements represent =b 1 standard deviation. Model results 
are lOyr climatological averages for SimCONT (red) and SimH20 
(blue). The dotted lines represent model uncertainty due to internal 
variability in the climate model computed as =b 1 standard deviation 
for n = 10 model run years, (b) isoprene mission temperate sites. 


water-stressed control run because drought-induced stom- 
atal closure drives higher canopy temperatures that stimulate 
the isoprene emission. At the North American sites (Harvard 
Forest and UMBS), the isoprene emission rate is a factor of 
3 lower in SimCONT than SimH20. At these sites, the mea- 
sured average diurnal cycle lies between the SimCONT and 
SimH20 climatologies. Interestingly, SimH20 authentically 
reproduces the GPP average diurnal cycle (Fig. 8c) at both 
sites (RMSE ~0.1 gCm -2 h _1 , R 2 >0.94). SimCONT GPP 
average diurnal cycle is too low by a factor of 2-3. It may 
be tempting to conclude that the water stress function in the 
control run is too strong at these locations, either because 
of errors in the climate model’s soil moisture simulation 
(too dry), and/or because of the model’s assignment of PFT- 
specific stress onset and wilting points. Discrepancies in the 
model’s PFT-specific photosynthesis and isoprene emission 
parameters provide another possible explanation. For exam- 
ple, using a local emission potential (7 S ) of 100 pgC g _1 h _1 
for the deciduous PFT, more than double the generic value 
applied in Yale-E2 (Table 1), isoprene emission simulated 
in LPJ-GUESS is within 5-10% of measured values at the 
UMBS site but only about half the magnitude of the mea- 
sured values at the Harvard Forest 1995 site (Ameth et al., 
2007). Increasing 7 S to 160 pgC g _1 h _1 for the Harvard For- 


Harvard Forest (June-July-August) 




0:00 2:00 4:00 6:00 8:00 10:00 12:00 14:00 16:00 18:00 20:00 22:00 0:00 


Fig. 9c. (a-c) Comparison between measured and modeled isoprene 
emission and GPP average diurnal cycles for the time periods in- 
dicated in Table 6. Measurements are black lines. Error bars on 
the measurements represent =b 1 standard deviation. Model results 
are lOyr climatological averages for SimCONT (red) and SimH20 
(blue). The dotted lines represent model uncertainty due to internal 
variability in the climate model computed as =b 1 standard deviation 
for n = 10 model run years, (c) GPP temperate sites. 


est site in LPJ-GUESS improved the agreement with obser- 
vations to within 20 %. In contrast, the JULES model used 
identical 7 S values for temperate broadleaf (45 pgCg -1 h -1 ) 
and tropical (24 jagCg -1 h _1 ) as those applied in this study, 
but overestimated the isoprene emission flux at 6 of the 9 
sites (UMBS, Harvard Forest, LV, Montmeyan, Manaus and 
Santarem) with the maximum discrepancy at La Verdiere 
(+236%). In JULES, model isoprene emissions in the mid- 
dle of the day and in the afternoons are generally higher than 
observed at all test sites. 

Figure 9a and b suggest that the model isoprene emission 
ramps up too quickly in the morning hours compared to mea- 
surements at both tropical and temperate sites. The observed 
rate of emission decrease in the early evening is better cap- 
tured by the model. The existence of a circadian rhythm in 
the PFT-specific fraction of electrons available for isoprene 
synthesis (e) provides one possible explanation, such that 
s is lower in the morning hours during the first onset of 
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Table 7. Statistical performance of model simulated isoprene average diurnal cycles at 9 sites in the tropical and temperate zones. 


Site 

Linear 

regression coefficient RMSE (mgC 
for 30 min data 

m -2 h -1 ) 

Maximum in average diurnal cycle 


SimCONT SimH20 

SimCONT 

SimH20 

Measurement 

SimCONT 

SimH20 

Costa Rica 

0.89 

0.90 

3.03 

3.14 

1.32 

7.31 

7.65 

Manaus 

0.93 

0.93 

1.78 

1.11 

6.41 

8.68 

6.94 

Santarem (wet) 

0.82 

0.81 

3.09 

3.20 

0.38 

6.93 

7.00 

Santarem (dry) 

0.92 

0.95 

1.44 

2.32 

3.41 

6.13 

8.17 

Bomeo/OP3 

0.96 

0.96 

2.57 

3.04 

1.65 

6.78 

7.77 

Harvard Forest (1995) 

0.67 

0.90 

1.33 

2.62 

7.40 

3.37 

10.23 

Harvard Forest (2007) 

0.74 

0.95 

0.81 

3.23 

5.21 

3.37 

10.23 

UMBS (5yrs) 

0.64 

0.89 

1.21 

2.44 

4.07 

1.71 

8.61 

La Verdiere 

0.82 

0.90 

2.43 

4.01 

4.01 

7.89 

11.24 

Montmeyan 

0.81 

0.84 

1.78 

2.48 

8.95 

9.40 

11.41 

Belgium (3 yr) 

0.87 

0.88 

4.22 

4.52 

0.77 

9.25 

9.77 


Table 8. Statistical performance of model simulated GPP average diurnal cycles at 2 temperate sites. 


Site Linear regression coefficient RMSE (gC m 2 h 1 ) Maximum in average diurnal cycle 

for 30 min data 



SimCONT 

SimH20 

SimCONT 

SimH20 

Measurement 

SimCONT 

SimH20 

Harvard Forest (2007) 

0.97 

0.97 

0.34 

0.09 

1.08 

0.45 

0.95 

UMBS (5yrs) 

0.79 

0.94 

0.36 

0.10 

1.02 

0.37 

0.99 


photosynthesis (Hewitt et al., 201 1). We will explore the use 
of a diurnally varying s in future site-level work. 

Based on linear regression (Table 7), the model demon- 
strates significant skill in capturing the observed variability 
in the 30 min average diurnal cycle across all sites (64-96 %). 
The model performance is best in the tropics (> 80 %) where 
there is little difference between SimCONT and SimH20. In 
the temperate zone, SimH20 demonstrates superior perfor- 
mance over SIMCONT (e.g. >90% versus 70-80%) in re- 
producing the diurnal average variability, especially at Har- 
vard Forest and UMBS. Essentially, the model is able to re- 
produce all of the observed variability in the average diur- 
nal cycle at Santarem (dry), Manaus and Borneo. SimH20 
reproduces the diurnal variability better than SimCONT at 
Harvard Forest (2007) and UMBS (all years). Despite the 
model’s overestimate of emission magnitude at Belgium (be- 
cause model PFT is deciduous in this location versus mixed 
temperate at the actual site), the model captures > 87 % of the 
measured variability in average diurnal cycle. 


5 Discussion and conclusions 

Isoprene emission is a quintessential quantity in chemistry- 
climate interactions. In this study, we have implemented a bi- 
ologically realistic photosynthesis-dependent isoprene emis- 
sion scheme into a global chemistry-climate model frame- 
work. The model provides a new tool that will allow us to im- 


prove understanding of the multiple interactions between at- 
mospheric chemistry, land ecosystem physiology and climate 
at regional and global scales that manifest through BVOC 
emissions. Such improvements are essential to provide ad- 
equate assessment of the climate and air quality benefits of 
mitigation strategies involving the short-lived climate forcers 
(Ameth et al., 2009). 

Overall, based on comparison with above canopy flux 
measurements, the model provides an authentic representa- 
tion of isoprene emission diurnal variability in tropical and 
temperate ecosystems but is less successful in reproducing 
the magnitude due to the assignment of PFT-specific iso- 
prene emission potentials, and possibly biases in the climate 
model’s internal meteorology. Current generation global cli- 
mate models are limited to 5-20 PFTs. To improve simula- 
tion of isoprene emission magnitude in global models, more 
measurements of the plant-to-plant species level variability 
in isoprene emission (the isoprene emission potential) are 
needed as well as a way to represent the greater species di- 
versity in models. In addition to accounting for leaf age ef- 
fects, improvements in simulating isoprene variability may 
be achieved by including climate-sensitive phenology, vari- 
able atmospheric surface CO 2 concentrations and the effects 
of ozone on plant physiology. 

In future work, we will use the model to assess the 
impacts of the photosynthesis-dependent isoprene emission 
on atmospheric chemical composition and the impacts of 
global change on the isoprene emission. We will explore the 
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Table Al. PFT-specific photosynthesis parameters used in Yale-E2 global carbon-chemistry-climate model. 


PFT 

Fcmax 

(jimolm -2 s _1 ) 

Nleaf 

(g[N] m -2 [leaf]) 

s* 

s wilt 

SLA 

(m 2 [leaf] kg - 1 [leaf]) 

Tundra 

21 

1.60 

0.50 

0.30 

2.25 

Grass 

26 

3.27 

0.30 

0.10 

11.7 

Shmb 

17 

2.38 

0.40 

0.22 

3.25 

Savanna 

30 

3.12 

0.65 

0.22 

5.1 

Deciduous 

30 

1.07 

0.50 

0.29 

8.3 

Evergreen 

42 

1.80 

0.50 

0.25 

5.9 

Rainforest 

54 

2.70 

0.60 

0.29 

9.9 

Crop 

30 

2.50 

0.45 

0.27 

6.36 


implementation of a leaf demography scheme (Caldararu et 
al., 2012) that will allow us to account for the effects of leaf 
age on the isoprene emission, which has implications for 
atmospheric chemistry, for instance, timing of the seasonal 
transition from VOC to NO x -limited ozone production. We 
will explore daily and seasonal average isoprene emission 
model performance using an off-line version of the vegeta- 
tion model that is driven by meteorology from the GMAO 
Modern Era-Retrospective Analysis (Rienecker et al., 201 1). 

We will extend the model to include a photosynthesis- 
based monoterpene algorithm (Schurgers et al., 2009). A ma- 
jor advantage of the vegetation biophysics in Yale-E2 is the 
high temporal resolution (30 min) that will allow additional 
treatment of monoterpene short-term storage in the meso- 
phyll not considered in any previous studies. The short-term 
storage may be important in simulating a correct monoter- 
pene emission diurnal cycle. 

This study has adopted an isoprene emission model that 
depends on the electron requirement for isoprene synthesis 
(Niinemets et al., 1999) necessitating an empirical param- 
eterization of the CO 2 -inhibition effect. Recently, a simple 
conceptual model of isoprene emission based on regulation 
of metabolism in the chloroplast has been proposed (Harri- 
son et al., 2013). This formulation does not require empiri- 
cal parameterization of the CO 2 -inhibition effect. The model 
needs further testing and development before application at 
the global scale but does hold significant promise of an even 
more faithful isoprene emission representation for the new 
generation Earth-system models. 
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